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NON-DETECTION OF THE PUTATIVE SUBSTELLAR COMPANION TO HD 149382 
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ABSTRACT 

It has been argued that a substellar companion may significantly influence the evolution of the 
progenitors of sdB stars. Recently, the bright sdB star HD 149382 has been claimed to host a sub- 
stellar (possibly planetary) companion with a period of 2.391 days. This has important implications 
for the evolution of the progenitors of sdB stars as well as the source of the UV-excess seen in 
elliptical galaxies. In order to verify this putative planet, we made 10 radial velocity measurements 
of HD 149382 over 17 days with the High Resolution Spectrograph at the Hobby-Eberly Tele- 
scope. Our data conclusively demonstrate that the putative substellar companion docs not exist, and 
they exclude the presence of almost any substellar companion with P < 28 days and M sin i > 1 Mjup. 

Subject headings: binaries: close - binaries: spectroscopic - planetary systems - stars: horizontal- 
branch - stars: individual (HD 149382) - subdwarfs 



1. INTRODUCTION 

The nature of the progenitors of subdwarf B (sdB) 
stars is not well understood. One theoretical formation 
method is through a common-envelop e (CE) phase of 
evolution with an orbital companion ()Han et all [20021 
|2003[) . As the progenitor star ascends the first red giant 
branch, interaction with a companion via a CE may cause 
it to lose large amounts of mass. When the progenitor 
star finally ignites core helium burning on the horizontal 
branch, it has a small radius, a high effective tempera- 
ture, and a low-mass hydrogen env elope that ca nnot sus- 
tain shell burning: it is a sdB star. iSokei] (|1998[ ) suggests 
that a substellar companion (a brown dwarf or planet) 
may be massive enough to survive this CE phase with- 
out merging with the host star. This implies that some 
sdBs may be created via the evolution of single stars with 
large planetary companions. In addition to their impor- 
tance as tracers of the many paths of stellar evolution, 
sdB stars have broader importance in that they may be 
the domi nant source of the UV-e xcess seen in elliptical 
galaxies (jBrown et al.lll997l[2008l) . 

HD 149382 is the brightest known sdB star, and thus 
provides an excellent opportunity to better understand 
the progenitors of sdB stars. It is an sdOB star with 
atmospheric parameters Teff = 35500 ± 500 K, \ogg = 
5.80 ± 0.05 (IGeier et al.l [2001 hereinafter G09), and a 
pattern of photospheric abund ances altered by diffusion 
processes (jBaschek et al.lll982l who found Toff = 35000 ± 
2000 K, \ogg = 5.5 ± 0.3). G09 claim to have detected 
a substellar companion to HD 149382 with a mass of 
8 — 23 Mjup using 16 radial velocity (RV) measurements 
collected from four different telescopes over a span of four 
years. G09 claim this is significant because it suggests 
that the progenitors of sdB stars can be stripped of their 
hydrogen envelopes by substell ar companion s during the 
CE phase, strongly supporting ISokeii (119981 ) . This is the 
first claim of a close substellar companion to an sdB star 
able to affect the evolution of its host. For brevity we will 
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use the word "planet" to refer to a substellar companion 
o f HD 149382, whethe r a true planet or a brown dwarf. 

[Jacobs et al.l (|2011l ) reported a non-detection of the 
claimed planet from their own RV measur ements of HD 
14938 2 using the HERMES spectrograph ([Raskin et al.l 
[20 lit ). Their upper limit on the RV variability of HD 
149382 was 0.8 kms^^, nearly three times lower than the 
amplitude reported by G09. Since a planet with the re- 
ported short orbital period would have important impli- 
cations, we independently sought verification of its exis- 
tence, as we describe here. 

2. DATA 

We observed HD 149382 on 10 nights between April 10- 
27, 2010 u sing the High Resolution Spectrograph (HRS; 
[Tulll998[ ) in its i? = 60, 00 configuratioii, at th e Hobby- 
Eberly Telescope (HET; iRamsev et al|[l99l . a 9.2m 
telescope stationed at the McDonald Observatory near 
Fort Davis, Te xas. We exploited the queue scheduling of 
this telescope ([Shetrone et al.l [20071) to efficiently make 
our 10 observations and acquire good phase coverage of 
the orbit. Typical single-exposure signal-to-noise ratios 
at blaze peak near 5500 A were 100 per 4-pixel resolu- 
tion element. The use of a single, fiber-fed spectrograph 
should reduce the potentially significant systematic ve- 
locity errors that G09 may have incurred combining mea- 
surements made from four different telescopes. 

Our observing procedure was to obtain two exposures 
of the star through the 2" fiber followed immediately by 
two exposures with the iodine absorption cell inserted in 
the beam. The cell served to measure night-to-night in- 
strumental offsets due to imperfect repeatability of the 
instrument configuration. We also obtained bracketing 
thorium-argon (ThAr) lamp exposures through the cali- 
bration fiber, and the standard ThAr exposures and flat 
lamp exposures through the iodine cell at the beginning 
or end of each night. On the first night, we obtained only 
one stellar exposure of each type, and on the final night 
wc obtained three. 

2.1. Reducing Spectra 
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The spectra were reduced using REDUCE, an ad- 
vanced image reduction package designed to reduce 
2-diniensional, raw echelle spectra into l-dimensional, 
wave length-calibrated spectra (jPiskunov fc Valentil 
I2002D . The reduction was modified slightly to accom- 
modate fiber-fed spectra, rather than slit spectra. Each 
observation contains information on two charge-coupled 
devices (CCDs), a "blue" CCD and a "red" CCD; we 
limited our analysis to the blue CCD, where we expect 
most velocity information to be located. We extracted 
47 orders of l-dimensional spectra from the blue CCD. 

When determining velocities, we stacked all stellar ex- 
posures from the same night into a single spectrum. The 
time between stellar exposures on the same night is al- 
ways less than 31 minutes; finding RV from a stacked 
spectrum for each night should not change our results, 
since any change in the RV of the star will be small over 
the "smoothing interval" of our exposures. Since we com- 
bined pure stellar exposures with stellar-plus-iodine-cell 
exposures, we were careful to analyze only those orders 
that contain no iodine lines. 

2.2. Visual Companion 

In addition to the putative substellar companion, 
HD 149382 fea t ures a known visual, stellar companion. 
I0stensen et al.l ()2005() detected this visual companion us- 
ing adaptive optics. The visu al companion is ab out 1" 
away (corresponding to 75 AU. [Jacobs et al.]|201lD while 
the HRS fiber diameter is 2", so we likely collected light 
from both stars. We therefore took care in our veloc- 
ity an alysis to select only lines expected from the pri- 
mary. lUlla fc Theilll (119981 ) detected an infrared excess 
in HD 149382, probab ly due to this visual companion, 
but IGeier et all (j2010D detect no spectral features from 
the companion in the optical. Since we also observed 
in the optical, and since any orbital motion between the 
stars must be too slow to be detectable over a 17 day 
span, we do not expect significant RV errors due to the 
companion spectrum. 

3. METHOD 
3.1. Line Selection 

To aid in selecting absorption lines that are unequiv- 
ocally from the primary star, and not from the faint 
companion, we generated a synthetic spectrum fr om a 
CasteUi & Kurucz modeQ ijCastelh fc Kurucz|[2003l) with 
the approximate parameters of HD 149382, as reported 
by G09. We started from a solar-abundance Kurucz at- 
mosphere model in LTE with T^s = 35000 K, log 5 = 5.0. 

We selected a hue list from the V ienna Atomic Line 
Database ('VALD: lKupka et al.|[2000t) spanning the wave- 
length range 415 — 500 nm, and correspon ding to the 
mode l selected. We used the SYNTH3 code (jKochukhovl 
I2007f ) to generate a stellar s pectrum, adopting abun- 
dances for elements derived bvlBaschek et al.l (119821 ) , and 
solar abundances of lAsplund et al.n2009[ ) for remaining 
elements. The resulting spectrum was broadened with a 
rotational profile of (usini — 2.4kms~^), and convolved 
with a Gaussian profile to match the R = 60, 000 HRS 
observation. 

The synthetic spectrum is not an exact match to the 
observed spectrum of HD 149382, as several lines are 
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Fig. 1. — An atlas illustrating the 15 lines used in this RV anal- 
ysis. Each plot presents a region of the spectrum centered on the 
indicated waveglength (listed in A) averaged over all observations. 
The span of the x-axis in each plot is given in Table [Tl The tick 
marks on the x-axes represent integer air wavelengths in A. 



present in the model but not the data (and vice versa). 
Wc note though, that such an exact match is not required 
for our purpose. The synthetic spectrum is of sufficient 
fidelity to unambiguously identify lines that belong to 
the primary star, and we chose these lines (i.e. present 
in the model and the data) to measure the RVs. 

We visually selected 15 regions containing moderately 
deep and narrow metal absorption features. Each of 
these 15 pixel groups contains one or more lines at wave- 
lengths predicted from the model. Our final selection 
of absorption features is given in Table [TJ and the aver- 
aged observed spectrum for each pixel group is shown in 
Figure [TJ 

3.2. Radial Velocity Determination 

A potential source of systematic error is the stabil- 
ity of the HRS wavelength solution. We checked our 
ability to determine precise wavelength solutions on dif- 
ferent nights by comparing the night-to-night shift in 
the wavelength solution to the shifts we measured in 
the iodine region of the stellar spectra. Since the io- 
dine exposures were made through the science fiber and 
have extremely high Doppler content, they provide a 
highly sensitive measurement of any instrumental drift 
affecting stellar velocity measurements. We compared 
the shifts measured in 192 63-pixel sections of the io- 
dine region (through a cross-correlation analysis in pixel 
space) to the shifts we measured in the same sections us- 
ing the nearest adjacent ThAr exposures (made through 
the calibration fiber). On each night, the median shift 
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TABLE 1 

Absorption Features Used for RV Measurement 



Min. A 


Max. A 


Median A 


Species 


Lines 


[A] (air) 


[A] (air) 


[A] (air) 




4241.33 




AO A 1 Q'7 


N II 


2 


4253.19 


4z54.iy 


4253.69 


T~^^^ TTT O TTT TT 

re III, b III, U II 


5 


42o4.5y 


4285.59 


4285. U9 


^ TTT O TTT TT 

Ca III, b III, vj II 


3 


4oUz.4z 


4303.51 


43U2.9 / 


/~1 „ TTT TT 

L^a ill, L> II 


3 


4332.27 


4333.38 


4332.83 


S III, O II, N III 


3 


4364.30 


4365.40 


4364.85 


S III 


1 


4399.16 


4400.25 


4399.71 


Ca III, Zn III 


2 


4446.67 


4447.65 


4447.16 


N II 


1 


4510.57 


4511.58 


4511.08 


N III, Ne II 


4 


4514.52 


4515.49 


4515.01 


N III, Ne II, C III 


3 


4529.96 


4531.08 


4530.52 


N II, N III 


2 


4606.81 


4607.88 


4607.35 


Ne II, N II 

N II, O III, Si IV 


2 


4630.11 


4631.28 


4630.70 


5 


4633.73 


4634.87 


4634.30 


N III, Ne II 


2 


4640.18 


4641.28 


4640.73 


N III 


1 



Note. — The absorption features predicted from our syn- 
thetic spectrum in our 15 pixel groups. Only the lines whose 
flux minima occur within the search area are included in the 
tabulation of species and number of lines. 

measured from all 192 sections differs by < 0.1 pixel 
between the ThAr exposures and the iodine region of 
the stellar spectra, indicating that the ThAr exposures 
track the instrumental shifts of the basement-mounted, 
temperature-controlled HRS as observed through the sci- 
ence fiber to better than ~ 100 ms~^. The stability of 
the wavelength solution (i.e., its night-to-night precision) 
is not a limiting factor in our RV precision. 

To determine a RV for each night's average spectrum, 
we rebinned our model spectrum into the wavelengths 
of each pixel group, using the wavelength solution deter- 
mined from the adjacent ThAr exposures. We then cal- 
culated the peak in the cross-correlation function (CCF) 
in pixel space between our model and our observation, 
and used the dispersion of the wavelength solution to 
calculate a velocity shift. The primary source of uncer- 
tainty in our RV measurements is the uncertainty in the 
location of the CCF peak due to the finite signal-to-noise 
ratio of the data. 

Each of the 15 pixel groups yielded its own set of RV 
measurements. We subtracted the bary centric motion 
of the telescope using routines provided by the Califor- 
nia Planet Survey, whic h are proven below 0.001 kms^^ 
(e.g. [Howard et al1l2010[ ). Over the timespan of our data, 
the change in barycentric correction amounts to almost 
7kms~^. Probably because of line blends, the CCFs 
for the different regions yield RVs that do not share a 
common apparent systemic velocity, although they in- 
depcntly detect the barycentric motion. Since we are 
primarily concerned with the acceleration of the star in 
its putative orbit rather than the systemic velocity of 
the star, we subtracted the mean RV of each pixel group 
from its time series. For each night, we then averaged 
the resulting differential RVs derived from each pixel 
group. We used the standard error of the mean of the 
relative RVs to determine the "internal" error on each 
night's measurement. Our RVs thus have an arbitrary 
zero point. We present our results after applying these 
corrections in Table [21 

4. RESULTS 



TABLE 2 
Radial Velocities 



Barycentric Julian Date 


Velocity 


[-2455000 days] 


[km/s] 


296.924 


-0.06 ±0.12 


298.915 


0.01 ±0.06 


302.906 


-0.13 ±0.07 


305.898 


-0.02 ±0.06 


307.878 


0.23 ±0.10 


308.884 


0.17±0.15 


309.896 


-0.24 ±0.15 


311.876 


0.10±0.10 


312.884 


-0.05 ±0.09 


313.900 


0.07±0.15 



The root-mean-square variation of our measured veloc- 
ities is 0.13 kms~^. The (reduced) xt value for the null 
hypothesis of no measured stellar acceleration is 1.61, 
suggesting that either our uncertainties are slight under- 
estimates or that the target is intrinsically variable at 
this level. 

To verify that these measurements are inconsistent 
with the claim of G09, we fitted our RVs to a circular Ke- 
pleria n model using the RVLIN code (jWright fc HowardI 
I2009D . with an initial period estimate equal to the value 
quoted by G09. Because we are testing a specific hy- 
pothesis, we limited the solutions to values within 5- 
sigma of the values quoted by G09, specifically P = 
{2.381, 2.401} days. We allowed all other orbital elements 
in the fit to float, including the phase, since the time 
elapsed since the G09 measurements is sufficiently long 
that the propagated phase is uncertain. 

Our best fit results are period P = 2.401 days and 
RV semi-amplitude K = 0.05 km s'^ (Fig. ^. This 
semi-amplitude is much smaller than the claim of 2.3 ± 
0.1 kms^^ by G09, and the difference lies well outside the 
uncertainties of the two measurements. We also fitted 
our data to a sinusoid with exactly the semi-amplitude 
and period reported by G09, but with unknown phase 
and zero pouit off set, using the IDL routine mpfitfun 
(jMarkwardtl |2009[) . Figure [3] demonstrates the ' discrep- 
ancy between our findings and those of G09. Broad- 
ening our search beyond a 5-sigma window around the 
putative period, we find that other best fit values for K 
for all similarly short periods are below our RV errors 
(Fig. m Table [2]) indicating that they are consistent with 
Okms~^. 

Because some pixel groups contain weaker lines than 
other groups, we repeat our analysis for subsets of the 
groups that omit the poorest features (e.g. omit groups 
9 and 12 in Table [1]). Our results do not change sig- 
nificantly for multiple subsets, indicating that our exact 
selection of lines is not important. Furthermore, as Fig- 
ure [3] indicates, since the absolute value of all the RVs in 
our pixel groups are well underneath the semi-amplitude 
found by G09, it is unlikely that any combination of pixel 
groups will provide a similar fit to that found by G09. 

We have no reason to suspect the existence of planets 
at other periods, and our experiment was designed only 
to confirm or disconfirm this particular claim. Nonethe- 
less, we have thoroughly searched our data for any pe- 
riodicity to determine the parameter space for planets 
that is excluded by our data. We restrict our search to 
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Fig. 2.— Best fitted phase-folded RV curve for HD 149382, limit- 
ing the period to a range near the previously reported value (gray 
area of Fig. \^ while allowing the semi-amplitude and phase to 
drift. The central value near Okms~^ does not represent the sys- 
temic velocity of the star, since we have subtracted the mean values 
of each line from their nightly values. Our best fit within 5-sigma 
of the period reported by G09 is P = 2.401 d and X = 0.05 km s'^, 
which is much less than the putative value oi K = 2.3±0.1 kms^^. 
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Fig. 3. — The gray boxes show our RV measurements for all 15 
absorption features on each night. The black asterisks show the 
mean of each night's measurements and the associated standard 
error of the mean. The thin line shows our best fit within a 5- 
sigma range of P near value reported by G09 (gray area of Fig. 

The thick line shows the sinusoid reported by G09. As in Fig 
[2] the central value near 0kms~^ does not represent the systemic 
velocity of the star. 

periods between 0.2-28 days, limited both by a semima- 
jor axis > 1 stellar radius as well as the longest period 
easily detectable with our data span. We further restrict 
our search to circular orbits. 

The observational cadence of the HET creates "blind 
spots" in our coverage near the harmonics of 1 sidereal 
day. We have excluded these regions from the following 
analysis (we define such a "blind spot" as any period 
for which we have no data points for a stretch of phase 
larger than 144 degrees, or 0.4 cycles, and for some higher 
harmonics where the data cluster tightly at two phases 
180 degrees apart). 

We determined the semi-amplitudes of the best fit si- 
nusoids for a well-sampled set of periods between 0.2- 
28 days. The largest semi-amplitudes we find are near 
0.17kms~^, which we find at a variety of periods across 
our search space. If we adopt the upper limit of the range 
of masses of HD 149382 reported by G09 of 0.53 Mq, 
then we can translate these semi-amplitudes into upper 
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Fig. 4. — Best fit semi-amplitude X as a function of period P. 
The thick, dotted, gray segment corresponds to a 5-sigma range 
centered on the period reported by G09, and includes the fit used 
in Fig. |3] At all points near the putative 2.391 ± 0.002 day period, 
our measured K is much less than the value of 2.3 ± 0.1kms~^ 
reported by G09. 




Period [day] 

Fig. 5. — Full period-minimum mass diagram for all periods 
probed in this work. "Blind spots" created by the cadence im- 
posed by the fixed elevation of the HET are shaded in gray; in 
these regions our phase coverage is too poor to put meaningful 
constraints on the presence of RV companions. Minimum masses 
correspond to the M sin i values for planets in circular orbits at the 
given period with the amplitude of the best-fit sinusoid to our data 
at each period. 

limits on the minimum companion masses (Msini) at 
each trial period. We find no best-fit companions with 
Msini in excess of 0.8Mjup, which makes the existence 
of any close-in giant planet with M sin i > I Mjup highly 
unlikely. Figure [5] illustrates our upper limits and blind 
spots. In principle, a more massive companion could ex- 
ist in a face-on orbit, or in one of our period blind spots 
(or perhaps in a highly eccentric orbit), but this remain- 
ing parameter space is extremely n arrow. 

Our results agree with those of I Jacobs et al.l (|2011| ). 
who also find no massive planets out to P = 50 days 
through an analysis of four He I lines, although our con- 
straints on the RV variability of HD 149382 are over four 
times more stringent. 

5. DISCUSSION 

There are several claims of substellar compani ons in 
orbit around hot subdwarf stars, summarized in iSchuhl 
((2OT0h . With the exception of HD 149382 and SDSS 
J08205-I-0008 (discussed below), these claims arc all 
based on the interpretation of residuals in an O — C 
diagram, in which observed ("O") moments of eclipse 
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or pulsation maximum are compared with a calculated 
("C") ephemeris. A quasi-sinusoidal pattern of residuals 
is interpreted in terms of a light-time effect across the 
orbit of the subdwarf. For the eclipsing systems, the in- 
ferred "pla net" is circumbinary ; for the case of pulsations 
(V391 Peg: ISilvotti et al.ll2007T ). the sdB star is otherwise 
thought to be single. Except for HD 149382 and CS 1246 
(see below) , all of the claimed orbital periods for substel- 
lar companions to sdB stars are long, in the range of 3 
to 16 years. Nearly always, a parabolic term is fitted to 
the O — C residuals along with the orbital sinusoid. The 
parabola represents either an evolving pulsation period 
or an evolving orbital period of the eclipsing binary, as 
angular momentum is transferred by tides from the or- 
bit to the binary companion, which undergoes "magnetic 
braking" . The relatively short observational timcspans, 
with the additional issue of correlations between model 
parameters describing the parabolic term and the light- 
time orbit, conspire to make these claims of substellar 
companions quite fragile — new observations tend to de- 
stroy the previous solution. (See [Parsons "eFaII[20l0l for 
examples of broken timing orbits in the case of eclips- 
ing binaries with white dwarf primaries.) Quasi-cyclical 
variation in O — C residuals for binaries may have causes 
other than orbital motion inv olving a third bo dy, such as 
the "Applegate mechanism" () Applegatd 1 19921 ) . in which 
mass quadrupole variations induced by stellar activity 
cycles alter the effective Kcplerian mass of the binary. 

The only timing-based {O — C) detection of a compan- 
ion to an sdB star that h as been confirmed b y an RV 
study is that of CS 1246 ()Barlow et al.ll2011al lH). Here 
the orbital period is 14.1 days, but the minimum mass 
of the companion is Msini w 0.13 Mq, above the sub- 
stellar limit. 

Some of the claimed substellar companions in long- 
period orbits may be corroborated with continuing ob- 
servations, and others will surely be found. It is unlikely, 
however, that these large orbits tell us anything directly 
about the form ation process of the sdB stars themselves. 
(jSchuhl (|2010f ) docs discuss the possibility that large or- 
bits may be indicators for the previous existence of inner 
planets, destroyed in a CE process.) It is also important 
to bear in mind that only a lower limit Msini to the 
companion mass M is found from these timing solutions. 

Short-period orbits of substellar companions around 
sdB stars are of more interest for elucidating the prior 
history of the sdB star. In addition to the idea that 
a planet may assist in re moving the hydrogen envelope 
from the proto-sdB stars ()SokeHll9^8l) , there is the possi- 
bility that "second-generation" planets may form in the 
debris left from the merger of two helium white dwarfs 
(the "helium planet" scenario oflsHvotti 20Q1). Note 
however, that since a merger between a helium white 
dwarf and a low-m ass hydrogen-burning s tar may also 
create an sdB star ([Clausen fc Wadell2011|) . a resulting 
second-generation planet need not be of exotic composi- 
tion. 

The formation of sdB stars through envelope stripping 
by planets is an intriguing hypothesis, but in light of the 
fact that the claimed planet HD 149382 b does not ex- 
ist, it currently has little observational support. To date 



the only such companion known is that in the cclips- 
in g sdB-Hbrown dwar f system SDSS J08205-h0008, found 
bv lGeier etlill (1201 1[ ). The orbital period is 0.096 days, 
and the estimated mass for the strongly irradiated com- 
panion, for reasonable estimates of the mass of the sdB 
star, is in the range 0.045 — 0.068 Mq, slightly below the 
hydrogen-burning limit. This pushes the lower mass limit 
for companions in short-period sdB systems to so mewhat 
lower masses than before fsee lGeier et al.|[201lL for dis- 
cussion of HW Vir systems), but hardly into the plan- 
etary realm. Given the potential importance of signif- 
icantly substellar companions to the understanding of 
the extremes of binary star evolutionary channels as well 
as the formation of second-generation planets, searches 
for such companions to sdB stars in short-period orbits 
should continue. 

Meanwhile, the putative plan et orbiting HD 1 49382 
has already been put to use bv lDe Marco et al.l ([20TT[ ) 
in revising the "a formalism" for CE binary evolution. 
Those authors make the CE efficiency parameter a a 
strongly decreasing function of mass ratio. This con- 
clusion relies strongly on the large value of a ~ 33 in- 
ferred for the putative HD 149382 "planetary" system. 
iDe Marco et al.l ([20111 ) recognized HD 149382 as an out- 
lier, and presented parametrized forms for a derived with 
and without including it in the sample of studied sys- 
tems. With our demonstration that this system does not 
exist (at least, as specified by G09), the weaker depen- 
dence of a on mass ratio is favored. 

6. SUMMARY 

Our findings conclusively refute the lGeier et al.l ([20091 ) 
claim of a substellar object orbiting the hot subdwarf st ar 
HD 149382. The orbital solution found bv lGeier etall is 
strongly inconsistent with our RV measurements. Our 
best fit results near the putative planet's best period are 
P = 2.401 days and K = 0.05 km s~^, which is consistent 
with Okms~^ given the uncertainties in our measure- 
ments. We can further confidently rule out planets with 
M sin i > 1 Mjup in circular orbits with periods < 28 
days, except for a few special cases arising from the ca- 
dence of our observations. 

We argue that the evolutionary origins of this hot sub- 
dwarf cannot involve an extant, close-in companion of 
giant planet mass or greater. 
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